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We consider a square-integrable semimartingale and investigate the con- 
vex order relations between its discrete, continuous and predictable quadratic 
variation. As the main results, we show that if the semimartingale has condi- 
tionally independent increments and symmetric jump measure, then its dis- 
crete realized variance dominates its quadratic variation in increasing convex 
order. The results have immediate applications to the pricing of options on 
realized variance. For a class of models including time-changed Levy models 
and Sato processes with symmetric jumps our results show that options on 
variance are typically underpriced, if quadratic variation is substituted for 
the discretely sampled realized variance. 

1 Introduction 

1.1 Variance Options 

Let X be a stochastic process, and let V he a partition of [0, T] with n + 1 division points 

= t[} < < • • • < t;^ = T. 



*The authors would like to thank Johannes Muhle-Karbe, Mark Podolskij, David Hobson and Walter 
Schachermayer for discussions and comments. 
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Then the realized variance of X over V is given by 

n 

RViX,V) = ^{Xt,-Xt,_y. (1.1) 

k=l 

In a financial market, let S be the price of a stock, a currency future, or another security, 
and let Xf = log{St/So) be its (normalized) logarithm. Then the realized variance (1.1) 

of X, with the increments — tk-i typically being business days, is a way to measure 
the volatility of the security price over the time interval [0, T]. Often realized variance 
is scaled by ^ such that it measures 'annualized variance' over a time interval of unit 
length. Since the uncertainty of future volatility is a risk factor to which all market 
participants arc exposed, many wish to hedge against it, while others are willing to take 
on volatility risk against a premium. To this end volatility derivatives have emerged, 
that allow market participants to take up positions in realized variance, see e.g. Carr 
and Lee [6] for an overview. Many of these derivatives pay, at the terminal time T, 
an amount f{jiRV{X,V)) to the holder, where / is the payoff function^. We call such 
derivatives options on variance. Typical choices for the payoff function are 

variance swap: f{x) = x — K volatility swap: f{x) = s/x — K, 
call option: f{x) = {x — K)''^ put option: f{x) = {K — x)'^ 

where K G M^o- For the variance swap, K may be chosen in such a way that today's 
fair value of the swap is zero; this choice of K is called the swap rate, and we denote it 
by s. The strike K for e.g. call options can then be chosen relative to the swap rate s 
by setting K = ks for some k G M^o- We refer to this choice as relative strike options. 
Note that the payoff functions listed above, with the exception of the volatility swap, 
are convex. In this article we will be concerned with variance options defined through a 
generic convex payoff function / : R^o — >■ ^■ 

1.2 Quadratic Variation and the Convex Order Conjecture(s) 

One of the cornerstones of the valuation theory of variance options is the replication 
argument of Neuberger [23] (see also e.g. Lee [19]), which states that a variance swap 
can be replicated by holding a static portfolio of co-maturing European options on S, 
while dynamically trading in the underlying S. The strength of Neuberger's replication 
argument lies in the fact that it is essentially model-free, up to the assumptions that (a) S 
follows a continuous martingale, (b) European options of all strikes and co-maturing with 
the variance swap are traded, and (c) that the realized variance (1.1) can be substituted 
by the quadratic variation [X, X]7- of X. Here, we are mainly interested in the last 
assumption, which is based on the fact that for any semimartingale X, and sequence of 
partitions ('P"')neN of [0,T], the realized variance RV{X,V"') converges in probability 



^ Other examples are so-called weighted variance swaps, which are volatility derivatives, but whose 
payoff cannot be written as f{j;RV{X,V)); see Lee [19]. 
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to the quadratic variation [X, X]t as the mesh of the partition tends to zero; i.e. 



RViX,^) ^[X,X]t, inprobabiUty (1.2) 

as mesh('P") — > 0. For options on variance with non-linear payoff, the static rephcation 
argument of Neuberger breaks down, but in specific cases Uke the Heston model, dynamic 
replication strategies involving European options or variance swaps can be derived (see 
Broadie and Jain [3]). In general, even if perfect replication is not possible, the arbitrage- 
free price at time zero of an option on variance is given by e~^'^K [/(^i?y(X, V))] where 
E [.] denotes an expectation under the risk-neutral pricing measure, and r the risk-free 
interest rate. Also for risk-neutral pricing and imperfect hedging, realized variance 
is frequently substituted by quadratic variation, since the latter is both conceptually 
and computationally easier to use, and eliminates the dependency on the nature of the 
partition V. See [5, 9, 17] for examples of this approach. All this raises questions on the 
qualities of the approximation 



E 



f[^RV{X,P) 



E 



that is (a) how precise is it? and (b) is there a systematic bias? While the precision 
of the approximation has been studied in the asymptotic limit n — t- oo and T — )■ 
(see Broadie and Jain [4], Scpp [27], Farkas and Drimus [13] rcsp. Keller- Ressel and 
Muhle-Karbe [18]), we are here interested in the existence of a systematic bias without 
asymptotics, i.e. for fixed and finite n and T. Numerical evidence given in Gatheral 
[14], Bxihler [5], Keller- Ressel and Muhle-Karbe [18] strongly supports such a bias and 
suggests that the price of a variance option with convex payoff, evaluated on discretely 
sampled variance, is higher than the price of an option with the same payoff, evaluated 
on quadratic variation. Prom this evidence we are led to the conjecture that 

E [/ {RV{X, P))]>E [f {[X, X]t)] , (1.3) 

for all convex functions /, or in other words that realized variance dominates quadratic 
variation in convex order. We call this statement the 'convex order conjecture' between 
discrete and continuous variance and write it in more concise form as RV{X,V) >cx 
[X, X]t.^ a shghtly weaker version is obtained if (1.3) is required to hold only for all 
convex, increasing functions /. We call this the 'increasing convex order conjecture'; it is 
equivalent to the realized variance dominating quadratic variation in increasing convex 
order, or more concisely RV{X, V) >icx [X, X]t. If E [RV{X, V)] = E [[X, X]t] then the 
convex order conjecture and the increasing convex order conjecture are equivalent, see 
Shaked and Shanthikumar [28, Thm. 4.A.35]. 

Similar questions can be asked about the relationship between the quadratic variation 
[X, X] and the predictable quadratic variation {X,X). For continuous semimartingales, 



^Note that the factor i will be absorbed into the function / from now on. 

^This conjecture has been initially conceived in discussions with Johannes Muhle-Karbe at the Bachelier 
World Congress 2010 in Toronto. 



3 



of course, {X, X) coincides with [X, X] . For discontinuous processes {X, X) is different 
from [X, X] , but sometimes more analytically tractable and has been used as a substitute 
for reaUzed variance, e.g. in Kallsen et al. [17] for exactly this reason. In the same 
article systematic underpricing of variance puts with predictable quadratic variation in 
comparison to quadratic variation has been observed, such that one may conjecture the 
relation 

E[/([X,X]t)]>E[/((X,XM] (1.4) 

for all convex functions /. Note that {X, X) is the predictable compensator of [X, X] such 
that under suitable integrability assumptions E[[X, = E[(X, and increasing 

convex order is equivalent to convex order. The main goal of this article is to prove the 
presented conjectures under certain assumptions on X and to outline the consequences 
for the pricing of options on variance. 

1.3 Strategy of the proofs and related work 

The convex (or increasing convex) order relation Y <cx Z {Y <icx Z) between two 
random variables is a statement about the probability laws of Y and Z, and thus not 
sensitive to the nature of the dependency between Y and Z. Nevertheless it is often a 
useful strategy to couple Y and i.e. to define them on a common probability space, 
which allows to use stronger and more effective tools, typically martingale arguments.^ 
For the proof of the convex and increasing convex order conjecture between discrete and 
continuous quadratic variation, we will couple RV{X,V) and [X, X]r by embedding 
them into a reverse martingale, as the initial and the limit law respectively. 

The idea for the construction of this reverse martingale comes from the literature on 
quadratic variation of Levy processes, or more generally processes with independent in- 
crements. Cogburn and Tucker [11] show that the realized variance of a process with 
independent increments over a sequence of nested partitions 'P" converges almost surely 
(and not just in probability) to the quadratic variation. The crucial step of their proof 
is to show that the sequence of realized variances over nested partitions forms a re- 
verse martingale when the underlying process has symmetric distribution. The almost 
sure convergence then follows from Doob's martingale convergence theorem and can be 
extended to arbitrary processes with independent increments by a symmetrization ar- 
gument. This technique can in fact be traced back to Levy [20] where a corresponding 
result for Brownian motion is shown. 

Finally let us remark, that convex order relations in the context of variance options 
have also been explored by Carr, Geman, Madan, and Yor [10], but with a very different 
objective. The authors observe that in an exponential Levy model St = S^e^^, where 
X is a Levy process without Gaussian component, the annualized quadratic variation 
j[X, X]^ is decreasing with t in convex order. Consequently, for each convex payoff / 
the term structure t i-)- E [/(j[X, X]t)] of variance options is decreasing in t. Also in 
the proof of Garr et al. [10] a reverse martingale is used: the process 1 1-)- i[X, X]^. 



"'See Lindvall [21] for the many uses of couplings in probability theory. 
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2 Definitions and Preliminaries 



We briefly introduce some definitions and properties that will be used in the results and 
proofs of Section 3. 

For a stochastic process X and a partition V of [0, T], the realized variance RV{X, V) 
of X over V has been defined in (1.1). In the case that RV{X,V) has finite expectation, 
we also define the centered realized variance by 

RViX, V) = RV{X, -P) - E [RV{X, V)] . 

Finally, to allow for certain generalizations of our results, the h-centered realized variance 
is defined, for any function h : M^o by 

RV''{X,V'') = RV{X,V'') - h(K[RV{X,V'')]y (2.1) 

Realized variance and centered realized variance can be regarded as the special cases 
h{x) = and h{x) = x. We use the same notation to define an /i-centered version of the 
quadratic variation [X, i.e. we define 

{X;X]'^ = [X,X]T-h(E[[X,X]T]), 

provided the expectation is finite. We will only be interested in /7-ccnterings where 
h : M^o — >• IK is Lipschitz continuous with Lipschitz constant at most 1; we denote the 
set of these functions by Lip^(IR^o)- 

A sequence (P")„eN of partitions is called nested, if all division points of are also 

division points of P". Note that by inserting intermediate partitions we can always 
consider ('P") as a subsequence of a sequence (P") of nested partitions, where each V"' 
has exactly n + 1 partition points. For the results in this article it will be sufficient to 
consider partition sequences of this type. Finally, the mesh of a partition V is defined as 
usual by mesh (P) = sup^,g|]^^ „|(ffe— tfc_i). As pointed out above i?y(X, P") — t- [X, XJt^ 
in probability, if mesh(P") — )■ for any semimartingale X. This results holds also for 
non-nested partitions, see e.g. Jacod and Shiryaev [15, Thm. 1.4.47]. 
Finally, let (Qn) be a decreasing sequence of a-algebras, and let X„ be a sequence of 
integrable random variables such that X„ G Gn- Then (X„) is called a reverse martingale 
with respect to Gn, if 

E[Xn-l\Gn]=Xn (2.2) 

for all n G N. Note that (Gn) is decreasing, and thus not a filtration. Similarly, (Xn) 
is called a reverse submartingale if (2.2) holds with '>' instead of '='. An important 
result is the following (cf. Loeve [22, 29.3.IV]): A reverse submartingale converges almost 
surely and also in to a limit X^o- Note that in contrast to the (forward) submartingale 
convergence theorem, no additional conditions on {Xn) are needed. If we define the tail 
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cr-algebra Qoo = r\k=i Qn-> then for any n G N the hmit can be represented as 

IE [Xn\ Qoo] = X^o- 

3 Results on discrete and continuous quadratic variation 

3.1 Reverse martingales from realized variance 

Let F,P) be a filtered probability space, where F satisfies the usual conditions 

of right-continuity and P-completeness. Let Y be an F-adapted cadlag process, and 
let be a P-complete tr-algebra such that % <Z Fq. We say that Y is a process with 
^-conditionally independent increments if for all < s < t the increment Yt — Ys is 
independent of Fg, conditionally on T-L. This definition includes time-changed Levy pro- 
cesses, and additive processes in the sense of Sato [26]. The H-conditional independence 
is equivalent to the assertion that 

E [f{Yt - Ys)Z\ ^] = E [f{Yt -Ys)\n]-^ [Z\ %] (3.1) 

for all bounded measurable / and bounded J^^-measurable random variables Z. See 
also Kallenberg [16, Chapter 6] for results on conditional independence and Jacod and 
Shiryaev [15] for processes with conditionally independent increments. 

Moreover, we say that a process Y has ^-conditionally symmetric increments, if 

E [f{Yt -Ys)\H]=^ [f{Y, - Yt)\ n] (3.2) 

for all t,s > and bounded measurable /. Alternatively we can use conditional chara- 
teristic functions to characterize conditional symmetry, i.e. a random variable X is 
^-conditionally symmetric if and only if 

E[exp{iuX)\n]=E[exp{-iuX)\'H], VmGM. (3.3) 

It will be helpful to know that a process with conditionally symmetric and independent 
increments^ is a martingale up to an integrability assumption. 

Lemma 3.1. Let X be an F-adapted process satisfying E [\Xt\] < oo for all t > 0. If X 
has H- conditionally symmetric and independent increments, then it is a martingale. 

Proof. Let t > s > and let .Z be a bounded J^g-measureable random variable. Using 
first conditional independence and then conditional symmetry of increments we obtain 

E [{Xt - Xs)Z] = E [E [(Xt - Xs)Z\ n]] = E [E [(Xt - X,)| H] • E [Z\ %]] = 
= E [E [(X, - Xt)\ H] ■ E [Z\ H]] = E [(X, - Xt)Z] 



®Here and in the following 'conditionally symmetric and independent' is a shorthand for 'conditionally 
symmetric and conditionally independent' 
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and conclude that E [{Xt — Xs)Z] = 0. Since Z was an arbitrary bounded J-'s-measurable 
random variable it follows that E \ {Xt — Xa) \ J-s\ = and hence that X is a martingale. 

□ 

The following Lemma establishes a reflection principle for processes whose increments 
are conditionally symmetric and independent. 

Lemma 3.2. Let Y he a process with H- conditionally symmetric and independent in- 
crements and let > 0. Then the process Y defined by 



— YsAf — {Yg — Ys/\t* ) 



(3.4) 



is equal in law to Y , conditionally on %. 

Moreover, RV{Y,V) = RV(Y,V) for all partitions V for which t^ is a partition point. 

Proof. Showing that Y is equal in law to Y, conditionally on 71, is equivalent to showing 
that for any sequence to < ti < ■ ■ ■ tn and bounded measurable / : M" — )■ R it holds that 



E [fiYt^ -Yt„...,Yt„- Yt^_,)\ = E f{Yt, J 



n 



(3.5) 



To show this equality it is sufficient to prove the equality of the ^-conditional charac- 
teristic functions of each side. By inserting intermediate points we may assume without 
loss of generality that = tm for some m G {0, . . . ,n}. Using the ^^-conditional inde- 
pendence (3.1) and the ?^-conditional symmetry of increments we obtain 



E 



E 



E 



exp yYl^'' {^^» ~ ^^^-Oy 

/ n 

exp i^ufe (Ft, -yt^_j 



n 



n 



n 



•E 



•E 



expf-i Uk{Yt^-Yt^_^) 

\ k=m+l / 

(n 
i Uk (it, - 

k=m+l 



n 



n 



E 



k=l 



n 



which proves (3.5). It remains to show that RV{Y,V) = RV(Y,V) for partitions V for 
which t* is a partition point. Denote the partition hy to <■■■< tn as before and assume 
that is the m-th partition point, i.e. = tm- Then 



k=l k=m+l 
n 

= J2i^t,-Yt,_y = RV{Y,V). 

k=l 



□ 
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Using this Lemma we show that we can construct a reverse martingale from the re- 
alized variances of a process with conditionally symmetric and independent increments, 
when the realized variances are taken over nested partitions. For a process with (un- 
conditionally) symmetric and independent increments this result has been shown by 
Cogburn and Tucker [11, Thm. 1]. Our result is a minor variation of the theorem of 
Cogburn and Tucker, but will be generalized in the next section. 

Theorem 3.3. Let X be a process with 71 -conditionally symmetric and independent in- 
crements and let ('P"')neN be a sequence of nested partitions of [0, T] such that 
K[RV{X,V'^)] < oo for all n G N. Then the realized variance of X over the parti- 
tions {V") is a reverse martingale, i.e. it satisfies 

E[RV{X,V''-^)\gn] = RViX,V''), (neN) (3.6) 

where 

g^ = nv a(^RV{X,V),RV{X,V+^),...y (3.7) 

Proof. Since (V"') is a nested sequence of partitions, we may assume without loss of 
generality that P""^ and differ by a single division point, which we denote by t*. 
Denote by a, b the two closest division points of V^~^, i.e. t* is inserted into the interval 
[a,b]. Let X be the process X reflected to the right of t* as in (3.4). By Lemma 3.2 X 
is equal in law to X, conditionally on H. Moreover, also by Lemma 3.2, RV{X,V'') = 
RV(X,V^) for ah k > n. This implies that 

gn = nya(^RV{X,V''),RV{X,V''+^),...^ = (3.8) 
= H V a(^RV{X, V''),RV{X, 7^"+^), . . . ) . 

To ease notation we abbreviate 7^„ := a(jlV{X, V),RV{X, 7^"+^), . . . ^ The next step 

is to calculate U ■.= E [RV{X, 7?"-i) - RV{X, V) \ gn] ■ By direct calculation we obtain 
that 

RV{X,V^-^) - RV{X,V^) = (Xh - Xt,){Xt* - Xa) (3.9) 
and conclude from the properties of conditional expectations that 

E [UHQ] = E [{Xb - Xt*){Xt* - Xa)HQ] , Vi7 G Q e 7^„. 

Using Lemma 3.2 and (3.8) we obtain 

E [(Xb - Xt*){Xt* - Xa)HQ] =K[E[{Xt,^ Xt*)(^t* - Xa)Q\ H] H] = 

= E [e [(lb - Xt*){Xt* - Xa)Q\ n\ i/] = -E [{Xb - Xt*){Xt* - Xa)HQ] , (3.10) 

and hence that E \UHQ] = for all if G ^, Q £ T^n- From = ^ V TZn we conclude 
that indeed 

E [Rv{x,r''~^) - RV{x,r'')\ g„] = c/ = o 
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showing the reverse martingale property (3.6). 



□ 



3.2 Coupling Realized Variance and Quadratic Variation 

To introduce the quadratic variation to the setting we now add the assumption that X is a 
semimartingale, such that RV{X,V^) [X,X]t in probabihty whenever mesh('P") — >■ 
0. Furthermore we require that X is a square-integrable semimartingale, i.e. that 
it is a special semimartingale with canonical decomposition X = Xq + N + A where 
A'^ is a square-integrable martingale and A has square-integrable total variation. This 
assumption implies in particular that supjg[o,T](^t ~ -'^o)^ is integrable, see Protter 
[24, Cli. V.2].^ Combining these assumptions with Theorem 3.3 immediately gives the 
following Corollary. 

Corollary 3.4. Suppose that X is a square-integrable martingale with H- conditionally 

symmetric and independent increments and that {V"')n£N is a nested sequence of par- 
titions such that mesh(7'") 0. Then RV{X,'P'^) [X,X]t a.s., E[[X,X]t] < oo 
and 

K[RV{X,V'')\goo] = [X,X]T (3.11) 
holds for any V and with g^='HV n„eN (^{RV{X, V"-)). 

Proof Since RV{X,V") < 2{n + 1) sup^gjQ 'j-j(Xt — Xq)'^ and X is square-integrable, 
we have that E[i?F(X,P")] < oo for all n G N. It follows from Theorem 3.3 that 
the sequence {RV{X,V^))nm is a reverse ^„-martingale. By the convergence theorem 
for reverse martingales this sequence converges almost surely and also in , and it 
remains to identify the limit. Since mesh('P") — )■ we have by [15, Theorem 1.4.47] that 
RV{X.,V^) converges in probability to the quadratic variation [X, We conclude 
that [X, X]t is in fact the almost sure limit of RV{X, V^) as n ^ oo, that E ]\X, X]^] < 
oo and that E [RV{X, V'')\ Qoo] = [X, X]t. □ 

In the case that we are only interested in linear expectations of realized variance and 
quadratic variation, the symmetry assumption on the increments of X can be dropped 
and we obtain the following. 

Corollary 3.5. Let X he a square-integrable semimartingale with %- conditionally in- 
dependent increments and let V : = to < ti < . . . < tN = T be a partition of [0, T]. 
Then 

n 

E[Rv{x,v)\n] = E[[x,x]T\n] + J2^i^t>'-^t^-i\'^]^- (^-i^) 

k=l 

In paHicular, E[RV{X,r)\'H] = E[[X, X]t\'H] for all paHitions P o/ [0, T] , if and only 
if the process X is a martingale on [0, T] . 



These assumptions are very natural in the context of variance options, where X is the log-price of 

a security S — S^e^ . Under the pricing measure is a local martingale and the semimartingale 
property of X follows automatically by [15, 1.4.57]. The square-integrability of X then merely ensures 
that the prices of variance options (in particular variance swaps) are finite. 
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Proof. Let Y be an ^-conditionally independent copy of X. Then Z = X — Y is 
a process with ^-conditionally symmetric and independent increments. Decomposing 
RV{X-Y,V) we get 

RV{X -Y,V) = RV{X, V) + RV{Y, V)- (3.13) 

n-1 
A;=0 

which can be bounded from above by 2{RV{X.,V) + RV{Y,V)). Since X is square- 
integrable, we have that E [RV{X,V)\ < oo and it follows that also E [RV{X - Y,V)] < 
oo. Thus Theorem 3.3 applies and we obtain that 

E [RV{X -Y,V)\ Goo] = [X-Y,X-Y]t = [X, X]t + [F, Y]t. (3.14) 

Taking expectations and combining (3.14) with (3.13) yields 

E [ [X, x]t| + E [ [y, y]t\ n] = ^ [RV{x, v)\n] + ^ [RV{y, v)\n]- 

n-l 
k=0 

Since X and Y have equal H-conditional distributions (3.12) follows. It is obvious that 
the sum in (3.12) vanishes if X is a martingale. Conversely, if the sum vanishes for 
any partition V, then E[Xt — X^j'^] = for any < s < t. By the ^-conditional 
independence of increments, we conclude that for any Z e Tg 

E [{Xt - Xs)Z] = E [E [Xt - Xs\ 'H]Z]=E [Xt - Xs\ U] ■ E [Z\ U] = 0. 

Hence E [X^ — X^l Jv,] = and X is a martingale. □ 

3.3 Results for semimartingales with symmetric jump measure 

Theorem 3.3 on processes with symmetric and conditionally independent increments is 
typically not suitable for the applications to variance options that we have in mind. The 
reason is that the risk-neutral log-price X of some asset can be a martingale only in 
degenerate cases, as already the asset price St = Soe-^* itself must be a martingale (at 
least locally). Assume for illustration that the risk-neutral security price St = S^e^^ is a 
true martingale and non-deterministic, then it follows immediately by Jensen's inequality 
that E [Xt\ < logE [St/Sa] = Xq and hence that X is not a martingale. Even if 5 is a 
(non-deterministic) strictly local martingale that is bounded away from zero, it follows 
by Fatou's lemma and Jensen's inequality that 

E [Xt] < hminf E [X^^rJ < logE [StArJSo] = Xo, 



Such a copy can be constructed using disintegration with respect to H. 
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for some localizing sequence (r„)„(=N, and hence that X is not a martingale. For this 
reason we want to weaken the symmetry assumptions on X. We introduce the following 
assumptions: 

Assumption 3.6. The process X is a square-integrable semimartingale starting at 
Xq = with the following properties: 

(a) X has ?^-conditionally independent increments; 

(b) X has no predictable times of discontinuity; 

(c) the predictable compensator of the jump measure of X, is symmetric, i.e. it 
satisfies ^{u, dt, dx) = u{u, dt, —dx) a.s. . 

Remark 3.7. 'No predictable times of discontinuity' means that /\Xr = a.s. for each 
predictable time r. This condition is also called quasi-left- continuity of X. It does not 
rule out discontinuities at inaccessible stopping times, like the jumps of Levy processes. 
For processes with unconditionally independent increments, 'no predictable times of 
discontinuity' can be replaced by 'no fixed times of discontinuity' (see [15, Cor. II. 5. 12]). 

Let us briefly summarize some of the consequences of these assumptions. First, since 
X is a square-integrable semimartingale, it is also a special semimartingale (cf. [15, 
II. 2. 27, II.2.28]), and we can choose h{x) = x as a truncation function, relative to which 
the semimartingale characteristics {B, C, u)h are defined. For this choice of truncation 
functions we have that B = A, i.e. the drift B is exactly the predictable finite vari- 
ation part of the canonical semimartingale decomposition. Second, a semimartingale 
has ^-conditionally independent increments, if and only if there exists a version of its 
characteristics {B, C, v) that is H-measurable. The increments are unconditionally in- 
dependent, if and only if % is trivial or equivalently if there exists a version of the 
characteristics [B^C^p) that is deterministic (cf. [15, Thm. II. 4. 15]). Third, the as- 
sumption that X has no predictable times of discontinuity implies that the version of 
the characteristics might be chosen such that they are continuous (cf. [15, Prop. II.2.9]). 

Theorem 3.8. Let X he a process satisfying Assumption 3.6, and let V he a partition 
of [0, T] . Then there exists a a-algebra Q^o C IF such that 

w.[Rv{x,r)\goo] = [x,x]T + RV{B,r). (3.15) 

// the increments of X are unconditionally independent, it also holds that 

> WMt^ (3.16) 
for each h G Lipi(M^o); O'f'd with equality for h{x) = x. 

Remark 3.9. In the proof below Qoo will be constructed explicitly as the tail cr-algebra 
of a sequence similar to (3.7). 



E 



RV {X, V) 
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Proof. Let ('P") be a nested sequence of partitions with mesh(P") — )■ 0, such that 
jjN _ -p £qj, gQjj^g jY. Let X = y + 5 be the canonical decomposition of X into local 
martingale part and predictable finite variation part. Since X is square-integrable Y is 
a square-integrable martingale and B has square-integrable total variation. In addition, 
define 

= ny (t{rv{y,v''),rv{y,v'"+^), . . .) (3.17) 



and let Qoo be the tail u-algebra HneN^" {Qn)n&i- Following [15, Thm. IL6.6], the 
conditional characteristic function of an increment oiY = X — B \s given by 



<l)t,s{u) ■= E 



^iu{Yt-Ys) 



n 



2 ft r . \ (3-18) 



= exp Y(<^t ~ '^*) + j[ / (^'"'' ~ ^ ~ ^«^)^('^' dr, dx)^ . 

The symmetry of v implies that 4>t,s{u) = ^t,s{—u), which by equation (3.3) implies 
the H-conditional symmetry of Yt — Y^. Hence K is a process with "H-conditionally 
symmetric and independent increments. We also have that E P")] < oo since Y 

is square-integrable. Thus Theorem 3.3 can be applied to Y and we obtain that 

^[RV{Y,V'')\goo] = [Y,Y]T, (3.19) 

for all n G N and in particular for n = TV and V"' = V. Since X is quasi-left-continuous, 
its predictable finite variation component B is even continuous, and it follows by [15, 
Prop. L4.49] that 

[X, X]t = [F, Y]t + [B, B]t = % Y]t. (3.20) 
Furthermore we can decompose the realized variance of X as 

n 

RV{X,V) = RV{Y,V) + RV{B,V) + ^ {Yt, - Yt,_,) ■ {Bt, - Bt,_,) . (3.21) 

fe=i 

We set ?7 := E [{Yt^ - Yt^_^) ■ [Bt^ - g„] and show that C/ = 0, similar to the 

proof of Theorem 3.3. We set 7^„ := a(^RV{Y,V'^),RV{Y,V'^+^), . . .) and using the 
reflection principle from Lemma 3.2 we obtain for all G H and Q G Tin that 

E [UHQ] = E [E [ {Yt^ J Q\ H] {Bt^ Bt^_^) H] = 

= E [E [ {Yt^_^ - Yt^) Q\ n] {Bt^ - i?t,_0 H] = -E [UHQ] 

and hence that indeed, f/ = 0. Taking ^co-conditional expectations, we have 

^[{Ytu-ytu-.)-{Bt,-Bt,_;)\g^]=Q, 

and thus 

¥.[RV{X,V)\goo]=^[RV{Y,V)\Goo]+RV{B,V). (3.22) 
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Combining (3.19), (3.20) and (3.22) yields (3.15). Suppose now that the increments of 
X are unconditionally independent. Then B is even deterministic (see the discussion 
after Assumption 3.6), and we derive from (3.15) that 



E 



RV^{X, V) I ^oo] = [X, X]t + RV{B, V)-f(E [[X, X]t] + RV{B, V)) . 



Applying the Lipschitz property of / (3.16) follows. 



3.4 Convex Order Relations 



□ 



By a simple application of Jensen's inequality, Theorem 3.8 can be translated into a 
convex order relation between the realized variance and the quadratic variation of X. 

Theorem 3.10. Let X he a process satisfying Assumption 3.6, and let V he a partition 
of [0, T] . Then the following holds: 

(a) For any increasing convex function f 

E [f[RV{X,V))] > E [/([X,X]t)] . (3.23) 

If X is a martingale, then (3.23) holds for any convex function f. 

(h) If the increments of X are unconditionally independent, then the h-centered realized 
variance satisfies 

(3.24) 



E 



f(Rv\x,r))]>W.[f([X,xfr) 



for any increasing convex function f and h G Lip]^(M^o)- If X is a martingale or 
h{x) = X, then (3.24) holds for any convex function f . 

Summing up the theorem, we can say that the 'increasing convex order conjecture' set 
forth in the introduction holds true at least for square-integrable semimartingales X with 
conditionally independent increments, symmetric jumps and without predictable times 
of discontinuity. If X is also a martingale, then even the 'convex order conjecture' holds 
true. However, as discussed at the beginning of the section, X cannot be a martingale 
if is assumed to be a non-deterministic martingale, such that only the increasing 
convex order conjecture holds in the relevant applications in finance. 

Proof. Applying an increasing convex function / to both sides of (3.15) and using 
Jensen's inequality, we have that 



E 



f{RV{X,V) 



>f[[X,X]T 



Taking (unconditional) expectations the convex order relation (3.23) follows for all in- 
creasing convex /. If X is a martingale, then by Corollary 3.5 E [RV{X, V)\ = E \\X, X\t\ 
and it follows from [28, Thm.4.A.35] that (3.23) holds even for any convex /. Equa- 
tion (3.24) is derived analogously from (3.16). In the case that X is a martingale or 
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h{x) = X 'increasing convex' can again be replaced by 'convex' since E RV {X, V) 



E 



[X,X], 



□ 



4 Results on predictable quadratic variation 

At this point we want to bring the predictable quadratic variation {X^X) of X into 
the picture. As discussed in the introduction wc conjecture the convex order relation 
[X, X]t >cx {X,X)rp for many processes of interest. Indeed the following holds. 

Theorem 4.1. Let X he a square-integrable semimartingale with Ti- conditionally in- 
dependent increments. Then the quadratic variation [X,X] dominates the predictable 
quadratic variation {X, X) in convex order, that is 



E 



for any convex function f. 



f([X,X]T)] >e[/((x,x)^) 



(4.1) 



Proof. Let {B, C, v) be the characteristics of the semimartingale X relative to the trun- 
cation function h{x) = x. Following Jacod and Shiryaev [15, Sec. II. 2. a] the quadratic 
variation of X is given by 



[X,X]t = C|+ f [ x'^fi{uj,dx,ds), 
Jo Jr 



where //(w, dx, ds) is the random measure associated to the jumps of X. The predictable 
quadratic variation in turn is given (cf. Jacod and Shiryaev [15, Prop. II. 2. 29]) by 



{X,X), = Cf + 



f [ x^v 

Jo Jr 



{lo, dx, ds). 



Furthermore, as X has ^-conditionally independent increments, the characteristics {B, C, h 
of X are "H-measurable. The process [X, X]^ — {X,X)^ is a local F-martingale, and due 
to the square-integrability assumption on X a true F-martingale. Together with the fact 
that H C To we conclude that 

E [ [X, X]t - {X, X)^\ ^] = E [E [ [X, X]t - {X, X)^\ T^] \n] = 0. 

Applying a convex function / to both sides of E[[X, = {X,X)j, and using 
Jensen's inequality we obtain 



E 



/([X,X]t)|h] >/({X,X)^). 



The desired result follows by taking expectations of this inequality. 



□ 
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5 Applications to variance options 



We now apply the results of Section 3 to options on variance. We distinguish between 
models with unconditionally independent increments, such as exponential Levy models 
and Sato processes, and models with conditionally independent increments, such as 
stochastic volatility models with jumps and time-changed Levy models. In the first class 
of models we obtain stronger results than in the latter. 

5.1 Models with independent increments 

5.1.1 Exponential Levy models with symmetric jumps 

In an exponential Levy model the log-price X is given by a Levy process with Levy triplet 
(6, a'^ ,v). Clearly X is a semimartingale with independent increments, and without fixed 
times of discontinuity. If the Levy measure v is symmetric and satisfies J x^i>[dx) < oo 
then Assumption 3.6 holds with trivial T-L and the results of Section 3 apply. We list 
some commonly used Levy models with symmetric jump measure; for definitions and 
notation we refer to Cont and Tankov [12, Ch. 4], or in case of the CGMY model to 
Carr, Geman, Madan, and Yor [7]: 

• the Black-Scholes model, 

• the Merton model with mean jump size, 

• the Kou model with symmetric tail decay and equal tail weights (A+ = A_ and 
P = l/2), 

• the Variance Gamma process with = 0, 

• the Normal Inverse Gaussian process with = 0, 

• the CGMY model with symmetric tail decay {G = M), 

• the generalized hyperbolic model with symmetric tail decay (/3 = 0). 

5.1.2 Symmetric Sato processes 

A Sato process is a stochastically continuous process that has independent increments 
and is self-similar in the sense that Xxt = A'^Xj for all X,t > Q and for some fixed 
exponent 7 > 0. Such processes have been studied systematically by Sato [25], and used 
for financial modeling for example in Carr, Geman, Madan, and Yor [8]. In Carr ct al. 
[10] the authors consider options on quadratic variation, when the log-price process X is 
a Sato process. Following [25], a Sato process is a semimartingale with the deterministic 
characteristics {r'h,t^'^ a"^ ,i't{dx)) where Vfidx) is of the form ft(-B) = f l^Bj{t'~'x)i'{dx) 
for some fixed Levy measure v. If v is symmetric with f x'^u{dx) < 00 then X satisfies 
Assumption 3.6 with trivial Ti. Referring to the terminology and notation of [8], this 
symmetry condition is satisfied for the VGSSD model with 6 = 0, the NIGSSD model 
with 6" = and the MXNRSSD model with 6 = 0. 
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5.1.3 Order relations between option prices 

For the models described in 5.1.1 and 5.1.2, Theorem 3.10 impUes the following order 
relations between the prices of options on variance: 

Fixed Strike Calls Consider a call on variance with a fixed strike, i.e. with a payoff 
function .x (x — K)^ where K > {). In this case Theorem 3.10 applies with 
/i = and /(x) = {x — K)'^ . Hence the price of a fixed-strike call on discretely 
observed realized variance is larger than the price of its continuous counterpart. 

Relative Strike Calls Let s be the swap rate, and consider a call with payoff function 
X {ks — x)'^ . The difference to the fixed strike case is that now the swap rate s 
also depends on whether discrete or continuous observations arc used for pricing. 
If < 1, that is if the call is in-the-money or at-the-money then Theorem 3.10 
applies with h = kx and f{x) = {x)~^. Hence, the price of a relative-strike in-the- 
money-call on discrete realized variance is larger than the price of its continuous 
counterpart. 

ATM Puts and Straddles Consider an at-the-money put with payoff x {s — x)'^ or 
an at-the-money straddle with payoff x i— t- |s — x|, where s is the swap rate. In 
these cases Theorem 3.8 applies with h{x) = x and f{x) = (— x)+ or f{x) = \x\. 
Hence, the price of an at-the-money put or an at-the-money straddle on discrete 
realized variance is larger than the price of its continuous counterpart. 

5.2 Models with conditionally independent increments 

5.2.1 Time-changed Levy processes with symmetric jumps 

Consider a Levy process L with Levy triplet {b,a,u) where the jump measure u is 
symmetric with J x'^v{dx) < oo, and a continuous, integrable and increasing process tj, 
independent of L. Letting Tt act as a stochastic clock, we may define the time-changed 
process Xf = L^^. The process X is a semimartingale with characteristics (6x4, art, I'Tt). 
Defining % = a{Tt,t > 0) and Tt = a{Xs,0 < s < t) we observe that X is a F- 
adapted square-integrable semimartingale with ^-conditionally independent increments 
and symmetric jumps, and hence satisfies Assumption 3.6. 

5.2.2 Stochastic volatility models with jumps 

Suppose that the risk-neutral log-price X has the semimartingale representation 

dXt = btdt + at dWt + xi^ {li{dt, dx) - ut{dx)dt) , Xq = (5.1) 

where is a F-Brownian motion and /x(dt, dx) a F-Poisson random measure with com- 
pensator i't{x)dt, and ht,(Tt,vt are H-measurable with % C Fq. Moreover, assume that 
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vt(dx) is symmetric and that 



dt < oo, 



such that X is square- integrable (cf. Keller- Ressel and Muhle-Karbe [18, p. 8f|). Then 
X satisfies Assumption 3.6 and the results of Section 3 apply. This setting includes the 
uncorrelated Heston model and extensions with symmetric jumps. 



5.2.3 Order relations between option prices 

For the models described in 5.2.1 and 5.2.2, Theorem 3.10 implies the following order 
relation between the prices of options on variance: 

Fixed Strike Calls Consider a call on variance with a fixed strike, i.e. with a payoff 
function x ^ {x — K)'^ where if > 0. In this case Theorem 3.10 applies with ^ = 
and f{x) = {x — K)'^ . Hence the price of a fixed-strike call on discretely observed 
realized variance is always bigger than the price of its continuous counterpart. 



6 Conclusions and counterexamples 

In this article we have shown that the increasing convex order conjecture for realized 
variance holds true in a large class of asset pricing models, namely under the assumption 
that the log-price X is a square-integrable semimartingale with conditionally indepen- 
dent increments, symmetric jumps and no predictable times of discontinuity. However, 
the numerical evidence given e.g. in Keller-Ressel and Muhle-Karbe [18] suggests that 
the conjecture may hold true under even more general conditions. In particular for Levy 
models with asymmetric jumps the conjecture seems to be valid, although this case is 
not covered by the assumptions of the article. It is not obvious - at least not to the 
authors - how the strategy of the proof can be extended to these cases, in particular 
how the symmetry condition on the jumps can be relaxed. On the other hand removing 
the assumption of conditionally independent increments easily leads to a violation of the 
convex order conjecture. A numerical counterexample for a stochastic volatility model 
with nonzero leverage parameter has recently appeared in Bernard and Cui [1]. Another 
counterexample can be given using a time-changed Brownian motion; it is based on 
suggestions by Walter Schachcrmaycr and David Hobson. 

Counterexample 6.1. Let B be an F-Brownian motion and define the stopping time 
r = inf {t : ^ [— 1, 1]}. The stopped process is a bounded martingale and can 
be closed at infinity by adding the terminal value B]^ = Br, which only takes the 
values ±1. Let / be an increasing, continuous and bijective function from [0, 1] to [0, oo] 
(such as t I— 7- tan(t^)). Define on the interval [0,1] the process X = ((-^^)/(t))(g[Q 
which is a continuous martingale w.r.t F' = (•^/(t))jg[o i]' Using the continuity of / 
we obtain [X, X]i = limt_).i[i?^, i?'^]^^^) = r. Choosing the partition V = {0,1} we 
get RV{X,V) = (±1)^ = 1 a.s. Using Jensen's inequality and the fact that E [r] = 1 
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we conclude that [X, X]i >cx RV{X,V). Since [X, X]i is truly random the opposite 
relation RV{X,V) >cx [-'^j-'^]! cannot hold true. This counterexample can easily be 
extended to other partitions of [0, 1] by introducing additional stopping times. 

The next counterexample shows that also the convex order relation [X, X]t >cx 
{X,X)t between predictable and ordinary quadratic variation shown in Theorem 4.1 
collapses if the assumption of conditionally independent increments is violated. 

Counterexample 6.2. Let be a random variable taking the values ±1 with probability 
^ each. In addition, let r be a uniformly distributed random time on [0, 1]. Define the 
cadlag process Zt := ^l{^<t}. This process is zero up to time r and then jumps to one 
of the values ±1. Denote by EI the natural filtration of Z; it is obvious that Z is a 
H-martingale. By direct calculation 

[Z, Z]t = ^^l{T<t} = l{T<t}- 

In particular [Z,Z]i = 1 almost surely. By Bielecki and Rutkowski [2] the 'hazard 
function' H corresponding to the random time r is given by H{t) = — log(l — t) and 
l{T<t} ~ H{t A r) is a martingale. We conclude that t ^ H{t A r) is the predictable 
compensator oit^ ^{r<t} hence that 

(Z,Z), = -log(l-(iAr)). 

In particular {Z,Z)-^ = — log(l — r), i.e. it is exponentially distributed with rate 1, and 
hence (Z, Z)^ >cx [Z, Z]i, while the reverse relation does not hold true. 
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